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INTRODUCTION 

Consider Bose-Einstein condensate (BEC) of alkali atoms in the F = 1 hyperfine state, 
elongated in x direction and confined in the transverse directions y,z by purely optical 
means. The dynamics of this assembly of atoms is described by a 3-component normal- 
ized spinor wave vector 4> = (<£i,4>o,<I>_i)^(jc,£) satisfying the multicomponent non- 
linear Schrodinger (MNLS) equation 01 HQ], which in dimensionless coordinates can 
be written down as: 

idt&±i + d^ ±1 + 2(|<D ±1 1 2 + |4> | 2 )4> ±1 = 0, 

id&o + <9 2 <J> + (2|4> +1 1 2 + |$o| 2 + 2|4>_i | 2 )cJ>o + 2<£ u <£ + i<D_i = 0, 

The second model which describes BEC with F = 2 hyperfine structure is a 5-component 
MNLS system: 

id t <&±2 + d?<f> ±2 + 2 ((4> t ,3>) - |4> ±2 | 2 ) <f>± 2 + 24>^ 2 4> + i4>_ 1 - ^ 2 ^o = 0, 
idt& ±1 + d*& ±1 +2 ((3> f ,3>) - |<*>±i| 2 ) &±i+2% l &- l & +l +% l Q$ = 0, 



(2) 



idr®o + d~<P () + 2 I (4> f ,4>) - ^|4>o| 2 ) <J> ( > + 2<f>,,<f> + - 2<f>, ) <f>+ 2 4>_ 2 - 0. 



where (3>^,3>) = L 2 = _ 2 l^/tl 2 - Both models allow Lax representations and therefore 
are integrable by the inverse scattering transform method D4LI2LI3I]. The Lax pairs have 
natural Lie algebraic structure which relates them to the symmetric spaces BD.I ~ 
SO(n + 2)/SO(n) x SO(2) with n = 3 and n = 5 respectively. From algebraic point 



of view this means that the potential Q(x,t) of L takes the form Q(x,t) = [J,X(x,t)] 
where X(x,t) is a generic element of the Lie algebra so (n + 2) and the constant element 
J = diag (1,0, ...,0— 1) is a specially chosen element of the Cartan subalgebra f) C Q. 
For more details see [Sill]. 

The present paper extends the results of H, IH] for the class of MNLS related to BD.I- 
type symmetric spaces, i.e. for any n. We briefly outline how the direct and inverse 
scattering problem for the Lax operator are reduced to a Riemann-Hilbert problem. Next 
we find that a simple change of variables can cast the above-mentioned MNLS into the 
Kulish-Sklyanin model (KSM) J6Q. We also apply Mikhailov reduction group method 
[0] and derive several new types of MNLS interactions. We derive also the constraints 
on the polarization vectors in the dressing factors that are imposed by the reductions. 
Finally we apply a proper modification (see [0, @]) of the Zakharov-Shabat dressing 
method [0, |9fl and derive the soliton solutions of the MNLS and of KSM in particular. 
Thus we obtain several new types of integrable vector MNLS and their soliton solutions. 

The majority of papers devoted to soliton equations analyze and solve the inverse 
scattering problem (ISP) for the relevant Lax operators using the typical (lowest di- 
mensional) representation of the corresponding Lie algebra. At the end of our paper we 
briefly compare the properties of the dressing factors in two of the fundamental repre- 
sentations of the Lie algebra so(2r). We also elucidate some additional issues considered 
in [0, OH such as the structure of the soliton solutions and the effect of additional Z2- 
reductions. 

MNLS EQUATIONS FOR BD.I SERIES OF SYMMETRIC SPACES 

MNLS equations for the BD.I. series of symmetric spaces (algebras of the type so(n + 2) 
and J dual to e\) have the Lax representation [L,M] = as follows 

L\jf(x,t,?i) = id x yr+ (Q(x,t) - XJ)\ff(x,t,X) = 0. (3) 
M\j/(x,t,X) = id t \j/+(Vo(x,t)+XVi(x : t)-X 2 J)\j/(x,t,X) =0, (4) 

Vi(x,t) = Q(x,t), V (xj)=mdJ ld ^ + ^[ a dJ 1 Q,Q(x,t)]. (5) 

where ad jX = [J,X] and ad J 1 is well defined on the image of ad j in g; 

/ q T \ 
Q=\ p s q , 7 = diag(l,0,...0,-l). (6) 
V P T *o J 

The H-component vectors q and p have the form 

q = (<?i,..., q n ) T , p = (pi,...,p n ) T , (7) 
while the matrix so = enters in the definition of so(n): 

Xeso{n), X + S^X T S^=0, SM = £ (-l)^ 1 ^^, (8) 



for n = 2r + 1 and 

sM = i(-iy+\E$ +1 _ s + El%_ StS ) (9) 

for n = 2r. By Egp above we mean n x n matrix whose matrix elements are (2?™ = 
5 s i8pj. With the definition of orthogonality used in (|8]) the Cartan generators = 

E^l — E^\_ k n+ i_ k are represented by diagonal matrices. 

The Lax pairs, related to the symmetric spaces SO(n + 2)/(SO(n) x 50(2)) have 
special algebraic properties. They are determined by choosing J = H\ to be dual to 
e\ E E' '. It allows one to introduce a grading in g, i.e. g = Qo © 0i so that: 

[X h X 2 }eQ , [X h Yi]eQ h [Y h Y 2 }eg , (10) 

for any choice of the elements Xi,X 2 E go and Yi,Y 2 E 0i. The grading splits the set 
of positive roots of so(n) into two subsets A + = Aq U where Aq contains all the 
positive roots of g which are orthogonal to e\, i.e. (a,ei) = 0; the roots in /3 E A^ 
satisfy (/3,ei) = 1. For more details see |5|]. 

In writing down the Lax pair © we made use of the typical n x n representation of 
so(n). The Lax pair can be considered in any representation of so(n), then the potential 
Q will take the form: 

Q(x,t)= £ (q a (x,t)E a +p a (x,t)E- a )- (H) 
aeA+ 

Next we introduce n-component 'vectors' formed by the Weyl generators of so(n + 2) 
corresponding to the roots in A^: 

E l = { E ±{ ei -e 2 )T ■ ■, E ±(e l -e r ), E ±e l ,E±( ei+er ),. ■ ■ , E ±(e l+ e 2 )), (12) 

for n = 2r + 1 and 

E f = ( E ±(ei-e 2 )i- ■ ■ , E ±{e l -e r ), E ±{e x +e r ), ■ ■ ■ , E ±{ ei +e 2 ))i (13) 

for n = 2r. Then the generic form of the potentials Q(x,t) related to these type of 
symmetric spaces can be written as sum of two "scalar" products 

Q(x,t) = (q(x,t)-E+) + (p{x,t) -E-). (14) 

In terms of these notations the generic MNLS type equations connected to BD.I. 
acquire the form 

iq t + qxx + 2(q 1 p)q- (q,s q)s p = 0, 
ipt ~ Pxx ~ 2(q, p)p + (p, s p)s q = 0, 

With the typical reduction p k = q* k it gives: 

iq t + qxx + 2(q\q)q-(q,s q)soq* =0, (16) 



If we put n = 3 and introduce the new variables 4>±i = q\ 3, 4>o = q% we recover 
equations ([I])- Likewise with n = 5 and <J>± 2 = ^1,5, ^±1 = q2,4, ^0 = #3 we find eq. 
©. The Hamiltonians for the MNLS equations (fT3T) are given by 

/oo / 1 \ 

^l(<9 x p ,d x q)-{p ,q) +-(p ,Sop)(q ,s q)\, (17) 



THE DIRECT AND THE INVERSE SCATTERING PROBLEM 



The fundamental analytic solution 

Herein we remind some basic features of the inverse scattering theory for the operator 
L ©, see 10,0]. There we have made use of the general theory developed in [2j], lTL IT2L 
ll3N and the references therein. The Jost solutions of L are defined by: 

lim (j)(x,t,X)e ajx = n, lim y(x,t, X)e ajx = 11 (18) 

and the scattering matrix T(X,t) = i// _1 ^)(jc, f, X). The special choice of / and the fact 
that the Jost solutions and the scattering matrix take values in the group SO(n + 2) we 
can use the following block- matrix structure of T(X,t) 

/ m+ -t>- T \ I mj B T c\ 

T{X,t)=\ b+ T 22 s B , T(X,t)=l -B+ f 22 s b- 
\ B +T s m\ J \ c| -b +T s m+ 

(19) 

where ^ ± (A,?) and B ± (X,t) are n-component vectors, T 22 (A) is n x n block matrix, 
and mf(X), and cf(X) are scalar functions. Such parametrization is compatible with the 
generalized Gauss decompositions of T(X) which read as follows: 

T(X,t) = Tj-D+Sj, T(X,t) = T+DJSJ, 

f? = AP&) i Sf = e±^), D± = diag ((mtf\mUmf) Tl )- 

The functions mf and n x n matrix- valued) functions rtij are are analytic for X G C±. 
We have introduced also the notations: 



+ ■ 



B~ _ B+ ^, b 

p =—, x =— , P =-, 

mj mj mj 

± m±(p ±r 5 p ± ) m^(f Tr 5 T T ) 

c — — 



1 2 2 

b- = ^-B-, b+ = ^B\ B+ = s °ti'* s °b + , B- = ^b-, 

where n + = — b + b~ ,T /(2mt), \x = — sqB + B~^ sq/ (2m~y). There are some 
additional relations which ensure that both T(X) and its inverse f(X) belong to the 
orthogonal group SO(n + 2) and that T(X)f(X) = 11. 



Important tools for reducing the ISP to a Riemann-Hilbert problem (RHP) are the 
fundamental analytic solution (FAS) x ± ( x J 1 X). We will introduce two pairs of FAS 
using the generalized Gauss decomposition of T(X,t), see [H, [3 14]: 



xH*,t,X) = (j>(x,t,X)Sf(t,X) = xir(x,t,X)T?(t,X)Df(X), 
X r,± (x,t,X) = <t>(x,t,X)Sf(t,X)Df(X) = yr(x,t,X)T?(t,X). 

More precisely, this construction ensures that ^(x, X) = X ± (x,X)e lXjx and E, f,± (x, X) = 
X f,± (x,X)e' Xjx are analytic functions of X for X G C±. If Q(x,t) is a solution of the 
MNLS eq. (fT5l) then the matrix elements of T(X) satisfy the linear evolution equations 

[ill 

± A 2 ^, X) = 0, i d -^- ± X 2 B ± (t,X) = 0, 

dt (21) 

am, dm-, 

i ~r L = °» '-7^ = °- 

a? at 
Thus the block-diagonal matrices D ± (X) can be considered as generating functionals of 
the integrals of motion. The fact that all (2r — l) 2 matrix elements of mt (X) for X G C± 
generate integrals of motion reflect the superintegrability of the model and are due to the 



for X G C± and that their zeroes and poles determine the discrete eigenvalues of L 



degeneracy of the dispersion law of (fT5l) . We remind that Dj (X ) allow analytic extension 



The Riemann-Hilbert Problem 

The FAS for real X are linearly related [0, 0] 

X + (x,t,X)= X -(x,t,X)G ,j(X,t), G ,j(X,t)=Sj{X,t)S+(X,t) 
x'< + (x,t,X)=x''-(x,t,X)G' 0J (X,t), G' 0J (X,t) = f J + (X,t)T J -(X,t). U ^ 

One can rewrite eq. (|22l) in an equivalent form for the FAS ^(xJ^X) = x ± {x,t, X)e lXjx 
and t;' >:iz (x,t,X) = x'' ± (x,t,X)e l ^ ,Jx which satisfy the equation: 



i^ + Q(x)^(x,X)-X[J,^(x,X)}=0, 
i^ + Q( x )^ '.± (x,X)-X [J, (*, X )] = 



(23) 



and the relations 



lim ^(x,t,X) = 11, lim^' ± (x,f,A) = ll. (24) 



Then these FAS satisfy the RHP's 

^ + (x,t,X) = ^(x,t,X)Gj(x,X,t), Gj(x,X,t) = e- aj ( x+ ^GJ(X,t)e iXj ( x+Xt \ 

^+(x,t,X) = ^-(x,t,X)Gj(x,X,t), Cf J {x i X i t) = e- ,u ^Cfj{X i t)^ J ^\ 

(25) 



Obviously the sewing function Gj(x,X,t) (resp. G'j(x,X,t)) is uniquely determined by 
the Gauss factors Sf(X,t) (resp. T^(X,t)). In addition Zakharov-Shabat's theorem [|8|] 
states that if sewing functions Gj(x,X,t) and G'j(x,X,t) depend on x and t in the way 
prescribed above ensures that the corresponding FAS satisfy the linear systems ((231) . 

Assume we have solved the RHP's above and know the FAS £ + (x,t,X). Then the 
corresponding potential of L is recovered by 

Q{x,t) = lim X (j-Z + (x,t,X)j£ + (x,t,X)) . (26) 



REDUCTIONS OF MNLS 

The reduction group proposed by Mikhailov [7] provides four classes of reductions 
which are automatically compatible with the Lax representation of the corresponding 
MNLS eq. 

The reduction group Gr is a finite group which preserves the Lax representation 
[L,M] = 0, i.e. it ensures that the reduction constraints are automatically compatible with 
the evolution. Gr must have two realizations: i) Gr C Autg and ii) Gr C ConfC, i.e. as 
conformal mappings of the complex A -plane. To each g k G Gr we relate a reduction 
condition for the Lax pair as follows QVll : 

C k (L(T k (X))) = r lk L(X), C k (M(T k (X))) = r] k M{k), (27) 

where C k G Aut $j and T k (X) G Conf C are the images of g k and r\ k = 1 or — 1 depending 
on the choice of C k . Since Gr is a finite group then for each g k there exist an integer N k 
such that g, k = 11 . In all the cases below N k = 2 and the reduction group is isomorphic 
to Z2. More specifically the automorphisms C k , k= 1 , . . . , 4 listed above lead to the four 
possible classes of reductions for the matrix- valued functions 

U{x,t,X) = Q(x,t)-XJ, V(x,t,X)=V (x,t)+XVi(x,t)-X 2 J, (28) 

of the Lax representation: 

1) CiCtftfaW)) = u(X), c^vHmx))) = v(X), 

2) C 2 (U T (k 2 (X))) = -U(X), C 2 {V T (k 2 (X))) = -V(X), 

3) C 3 (U*(k 1 (X))) = -U(X), C 3 (V*(ki(A))) = -V(X), 

4) C,(U(k 2 (X))) = U(X), C 4 (V(k 2 (X))) = V(X), 



(29) 



In what follows we will examine the typical reductions of MNLS eqs. of the class 1) 
obtained by specifying K\(X) = X* and C\ to be a Z2-automorphism of q such that 
C\ (J) = J. Below we list several choices for C\ leading to inequivalent reductions: 



a) Ct = 11 , p(x) = q*(x), b) C\=K U p(x) = K 01 q* (x) , 

c) Ci=5 e2 , p(x) = K 02 q*(x), d) Ci=5 e2 5 e3 , p(x) = K 03 q* (x), 



(30) 



where 

Kj = block-diag (1,^,1), K 0l = diag (ei , . . . , e r _i , 1 , e r _i , . . . , £1) , (31) 



and £j = ±1. The matrices K02 and ^03 corresponding to the Weyl reflections S e2 , S e2 S e3 
etc. are not diagonal; they have dimension n x n and for n = 3, 4 and 5 are given by: 



3, Kq2 



n = 4, K 02 
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(32) 



5, Kq2 



Each of the above reductions impose constraints on the FAS, on the scattering matrix 
T(X) and on its Gauss factors Sj(X), Tj{X) and Dj(X). These have the form: 

{^{x,X*yf = Kj l $(x,X)Kj, { W {xX))^ = Kj l f{x,X)Kj, 
{ X + {x,X*)y = Kj l X-{^X)Kj [T{X*))^ = KjH{X)K h 
(5 + (A*)) f =Kj l S-{X)Kj (T+(X*))^ =Kj l f-{X)Kj 



(33) 



where the matrices Kj are specific for each choice of the automorphism C\, see eq. yT 
In particular, from the last line of (|33l) and (|3TI) we get: 



(m+(X*)r=m^(X), (34) 
and consequently, if m^(X) has zeroes at the points Xf, then mj~(A) has zeroes at: 



i. 



,N. 



(35) 



Below we will write down the effects of these reductions on the corresponding 
Hamiltonians. For the typical reduction p = q* we get: 



#MNLS 



dx{ ((d x q,d x q*)-(q,q*) 2 + l-(q,Soq)(q*,Soq 



(36) 



#MNLS = j_J x [id x qKjd x q*) - (q.Kjq*) 2 + ^(q,s q)(q*,s q*) S J , (37) 

The Hamiltonian #mn LS w ^ tn ^TT) has indefinite kinetic term. As a consequence the 
corresponding MNLS has singular soli ton solutions which 'blow-up' in finite time. 



The above Hamiltonians, after the change of variables can be written in more 'aes- 
thetic' form. Indeed, for odd n = 2r — 1 we can put: 



q2k-l,2r-2k+l 



V2&-1 ±n>2r- 2k+l 



Q2k,2r-2k 



W2k =F V2r-2k 



V2 ' V2 

with k= 1,2, . .. ,r — 1 and co,r = e^ -1 ) 7 "'/ 2 ; for n = 2r we put: 

V2fe-1 ±iV2r-2jfc+2 



C0,rV,-; 

(38) 



<?2fc-l,2r-2fc+2 ~ 



W2kTV2r-2k+l 
q2k,2r-2k+l ~ ^ w^J 



' ^*-T 1 ^ 

with & = 1,2, ... , r. 

Inserting the above changes of variables into the Hamiltonian (l36l) we get 



#KS 





2 










71 



















(40) 



which is the Hamiltonian of the n-component Kulish-Sklyanin model (KSM) [6]. Thus 
we have demonstrated that the Lax pairs ©, © can be used also for integrating 
the MNLS (l40l) . In their original paper Q6Q Kulish and Sklyanin have used Lax pair 
whose potential is an element of a Clifford algebra. Later Sokolov and Svinolupov 111 511 
discovered another class of Lax pairs for these models whose potentials take values in 
Jordan algebras. The above Lax pairs allowed to prove integrability of the KSM but 
were not convenient for solving the inverse scattering problem and constructing exact 
solutions. Another important property of these models is that they possess both classical 
0] and quantum 7?-matrices 

Another way to obtain KSM is to apply the reduction of type 4) with 
Kq = block-diag (l,£So, 1). where £ = ±1. For odd values of n = 2r — 1 this reduction 
means that: 

fc=l,...,r, (41) 
w k , k=l,...,r. (42) 



Qk = (-l) k+l £(l2r-k = Wk, 



while for n = 2r one gets: 

q k = (-l) k+l £q 2 r-k+l 



This reduction leads to r-component KSM. 

Let us write down the Hamiltonians for the different reductions. Below for conve- 
nience we will split //mnls into kinetic and interaction terms: //mnls = —H^. 



Reduction b): 



^kin 



H 



int 



J dx j £ £j ( I d x qj | 2 + | d x q 2r -j | 2 ) + | Mr 1 2 J 

/oo [ r— 1 

dxl Y,£j(\<lj\ 2 +\<l2r-j\ 2 ) + \qr\ 2 ) 2 



(43) 



I(-l)^ +1 2 W r-y+(-l)V 



One can construct other reductions, e.g. ones of type c) with reduction matrix Kj. 
Then 

#MNLS = f~JxUd x ^Kjd x q) - (fKjq) 2 + X - \(fs Q q)\ 2 \ , (44) 

Characteristic feature of the reductions involving Weyl group elements is that they 
lead to 'non-diagonal' form of the kinetic terms Biol . Making simple change of variables 
diagonalizing Kj we can recover the diagonal form of the kinetic terms but unfortunately 
we can not make it positive definite. This is related to the fact that Kj = 1 and so has as 
eigenvalues both +1 and —1 with certain multiplicities. 

Let us give also an important example of class 2) reductions (|28T) . The constraints that 
these class of reductions impose on the FAS and on the scattering matrix T(X) and on 
its Gauss factors Sf(X), Tj°(X) and Df (X) take the form: 

(Z + {x,l)) T = K'f 1 r {T{X)) T =K'fH{X)K' p 

(45) 

(S ± (X)) T = K': S^{X)K'j (r ± (A)) r = K'f x f ^(^K'j 

and (D ± (X)) T = K'-~ D^(X)Kj. The explicit form of the matrices Kj is determined 
by the particular realization of the automorphism C 2 . Choosing n = 3 and C2 = S ei we 
obtain the constraint qi= q-} and the reduced Hamiltonian takes the form: 

//=^Jx||^v 1 | 2 + |^v 2 | 2 -(|vi| 2 + |v 2 | 2 ) 2 + ^|v?-vi| 2 }, (46) 

where we have put q\=qi = 4|Vi and g 2 = v 2 - This model also been derived as relevant 
for F = 1 BEC |fl7h. 



DRESSING METHOD AND SOLITON SOLUTIONS 

The dressing Zakharov-Shabat method for constructing soliton solutions of MNLS 
has been modified in yO for the BD.I-type symmetric spaces. There we also analyzed 
the different types of soliton solutions. Below we briefly discuss the properties of the 
generic one-soliton solutions 

It is obtained by dressing the regular FAS %q(x, X) of the RHP (|25l) . Using them we 
construct the singular solutions Xq{x,X) of the RHP 

X ± (x,X) =u(x,X)Xq(x,X)u-, x L± ix,X) = u(x,X)Xq(x,X)u + , 
u{x,X) = ll + (ci(A)-l)/ 5 i(x) + (c7 1 (A)-l)AW, m±= lim u{x,X). (47) 

For the above choice of J it is enough to consider rank 1 projectors Pi (x, t) and P\ (jc, t) = 
SqP[Sq. Together with the constraint P\P\ = 0, the last condition ensures that u(x,t) G 
SO(n + 2). It remains to only to give the explicit form of P\(x,t). Generically it is 
determined by two polarization vectors |«o,l) an d (mo i|, and the initial regular solutions: 



fiM = ^^## ; P lM = l mM){n & , (48) 



\ni{x,t))=Xo{x,t,ki)\no,i), (*»i(*,OI = ( m o,i|2 ( x ^^i )> (49) 

=^o"(*,^V)l^oj)> <mi(x,OI = (noA\Xo(x,t,Xi). 
The one soliton solution is parametrized by the two eigenvalues Aj and by the polar- 
ization vectors |«o.i) an d (moi\- The latter after renormalization have n — l independent 
components each: 

( ^ \ f r- r- r-\ 

l*o,i) = v ,i/v^o , ( m o,i| = V B o,Mo,i/v 5 o,l/V 5 o ), 

V i/v^ / v ) 

where Ao = jC^oVoVo.i) and Bo = \{fil i*oPo,l )• The constraint Pi Pi = means that 
the vectors /lo.i and Vo.i must satisfy jUq i^oVo.i = 0- Therefore the one-soliton solution 
can be viewed as a dynamical system with In — I degrees of freedom. After some 
simplifications it takes the form: 

q k (x,t) = — — e [cos(^)cosh(z 0fc ) + /sin(^)sinh(z /t)] , (50) 

A = 2cosh(2z ) + e, C = " , , z k = x + w l t-8 . k /nu 

\ A o\ 

Zo = Vi(x-Uit) + fy), Zok=VlX+fy) jk , Zok= Vi(x-Vif)+|o } jfc, (51) 

e 1 -| I 1 1 e 1 , l v 0,2r-& g 1 , 

§o = -lnAo, go,* = ^ 0fe =- m 

2 2 |Ao||Vo,jt| 2 

5o,fe = («o,2r-yt + ao,fe -«o- rc*)/2, 8 0jk = (a 2r- k - ao,k + cc + nk)/2, 
where Oq = arg Aq and ao,A: = arg Vo,fc. 

Each of the reductions of the type (l29l) imposes constraints not only on Aj + = (Aj - )*, 
but also on the polarization vectors: 

00,1=^0,^0,1*, vJi^oj5oV ,i = 0. (52) 

As a result, after the reduction the number of independent parameters of the soliton 
solution becomes n—l. The velocities u\ and w\ are given by u\ = —1\X\ and w\ = 
(Vi 2 -M?)M- 

Special attention deserves the fact that generically all zo,k are different and as a result 
each component q k (x,t) has its center of mass shifted with respect to the others. 

Let us now consider a Z2 x Z2 reduction by applying simultaneously two reduction: 
the first is the typical one and the second is the class 2) reductions as for the model (1461) . 
The first reduction imposes the relation (|52l) between the two polarization vectors |«o,i) 
and («o,i I- The second reduction imposes constraint on the vector |«o,i)> namely: 

|*0,l) =K^\n 0>1 ). 

In particular, for n = 3 and C2 = S ej the vector |no,i) has 3 components and 
Kqj = diag (1, — 1, 1). Thus only two independent complex coefficients are enough 
to parametrize the corresponding polarization vector, and the corresponding soliton can 
be viewed as dynamical system with three degrees of freedom. 







|Aol 





DISCUSSION AND CONCLUSIONS 



One of the important consequences of the FAS is that with their help one can construct 



the kernel of the resolvent of L (see |12l . LL80) and prove the completeness relation 



for its eigenfunctions. From these expressions it becomes obvious that the resolvent 
develops poles at all points Xj^ G C± for which mf(Xj^) = 0. Combining this fact with 
the equivalence between the solutions of the RHP and the FAS of the Lax operator we 
conclude that the singularities of the RHP correspond to the discrete eigenvalues of L. 

Quite often the general analysis of the MNLS (OQ) is followed by simplifications which 
often reduce the MNLS to a single-component NLS. One way do to this was mentioned 
above: it is to impose the reduction 4> + i = 4>_i. Another less obvious way to this is 
to impose this reduction on the initial conditions. Indeed, one can show that impos- 
ing <3> + i(x, f = 0) = <t>-i(x,t = 0) ensures that & + \(x,t) = <J>-i(*,f) for all t > 0. At 
the same time there is a substantial difference between the solitons of the scalar NLS 
or Manakov model and the solitons of MNLS ([I])- Unlike the solitons of the Man- 
akov model, all three components of the one-soliton solution of (QQ) have different x- 
dependence; generically each component has different 'center of mass' position. There- 
fore, if one wants to demonstrate new nontrivial aspects of soliton dynamics one should 
use generic initial values for <fr±\(x,t = 0)and 4>_i(*,? = 0). 

Another still open problem is the interrelation between the solutions of the direct 
and inverse scattering problem for L, considered in different irreducible representations 
(IRREP) of the corresponding Lie algebra q. From the point of view of the relevant 
NLEE, their Lax representations have purely algebraic nature and therefore, the form of 
the NLEE does not depend on the choice of the IRREP of g. 

From the point of view of the spectral theory, the different IRREP have different 
dimensions; therefore changing the IRREP we change the order of the corresponding 
operator. Since we are dealing with simple Lie algebras whose IRREP are well known 
Q5D . In particular, it is well known that the finite dimensional representations can be real- 
ized as invariant subspaces of the tensor products of the typical one. Let us assume that 
we are able to construct the FAS and the relevant RHP and dressing factors in the typical 
representation. Obviously, taking the tensors products of the FAS their analyticity prop- 
erties will persist and we will get the corresponding FAS and RHP in the corresponding 
IRREP. However nontrivial things may take place when one considers the multiplicities 
of the corresponding discrete eigenvalues. 

As an example I will just mention that the dressing factor can be evaluated also for the 



other fundamental representations of g II 1911 . If in the typical representation of $j ~ so(2r 



u(x,X) is given by (1471) then in the spinor representation it will take the form [20]: 



1 



u(x,X) = y/c l (h)ni(x,t) + —==7C l (x,t), Tt\{x,t) = s Q ni(x,t)s 

c\{X) 



and the projectors satisfy 7t\(x,t)7i\(x,t) = and 7li(x,t) +7Z\(x,t) = 11. Note the sub- 
stantial change in the A -dependence of u(x,X), as well as the fact that now instead of 
having rank one projectors Pi (x, t) we get projectors %\ (x, t) and Tt\ (x, t) of rank r. 
We will discuss these problems in more details elsewhere. 



ACKNOWLEDGMENTS 



It is my pleasure to thank the organizers of the AMITANS conference for kind invitation. 
I am grateful to professor N. Kostov and to an anonymous referee for useful suggestions. 

REFERENCES 

1. J. Ieda, T. Miyakawa, and M. Wadati. Exact analysis of soliton dymamics in spinor Bose-Einstein 
condesates. Phys. Rev Lett. 93, 194102 (2004). 

2. V. S. Gerdjikov, N. A. Kostov, T. I. Valchev. Solutions of multi-component NLS models and Spinor 
Bose-Einstein condensates Physica D 238, 1306-1310 (2009); arXiv:0802.4398 [nlin.SI]. 

3. V. S. Gerdjikov, D. J. Kaup, N. A. Kostov, T. I. Valchev. On classification of soliton solutions of 
multicomponent nonlinear evolution equations. J. Phys. A: Math. Theor. 41 315213 (2008) (36pp). 

4. A. P. Fordy, and P. P. Kulish. Nonlinear Schrodinger equations and simple Lie algebras. Commun. 
Math. Phys. 89, 427^143 (1983). 

5. Helgasson S. Differential geometry, Lie groups and symmetric spaces, Academic Press, (1978). 

6. P. P. Kulish, E. K. Sklyanin. O (N) -invariant nonlinear Schrodinger equation - a new completely 
integrable system. Phys. Lett. 84A, 349-352 (1981). 

7. Mikhailov A V. The Reduction Problem and the Inverse Scattering Problem. Physica D, 3D, no. 1/2, 
73-117 (1981). 

8. V. E. Zakharov, A. B. Shabat. A scheme for integrating the nonlinear equations of mathematical 
physics by the method of the inverse scattering problem. I. Functional Analysis and Its Applications, 
8,226-235 (1974). 

V. E. Zakharov, A. B. Shabat. Integration of nonlinear equations of mathematical physics by the 
method of inverse scattering. II. Functional Analysis and Its Applications, 13, 166-174 (1979). 

9. V. E. Zakharov, and A. V. Mikhailov. On The Integrability of Classical Spinor Models in Two- 
dimensional Space-time Comm. Math. Phys. 74, 21^-0 (1980). 

10. Nikolay Kostov, Vladimir Gerdjikov. Reductions of multicomponent mKdV equations on symmetric 
spaces of Dill-type. SIGMA 4 (2008), paper 029, 30 pages; arXiv:0803.1651 . 

11. Zakharov V E., Manakov S V, Novikov S P., Pitaevskii L I. Theory of solitons. The inverse scattering 
method, Plenum, N.Y. (1984). 

12. VS. Gerdjikov. On the spectral theory of the integro-differential operator A, generating nonlinear 
evolution equations. Lett. Math. Phys. 6, n. 6, 315-324, (1982). 

13. V. S. Gerdjikov. Generalised Fourier transforms for the soliton equations. Gauge covariant formula- 
tion. Inverse Problems 2, no. 1, 51-74, (1986). 

14. V. S. Gerdjikov. The Generalized Zakharov-Shabat System and the Soliton Perturbations. Theor. 
Math. Phys. 99, No. 2, 292-299 (1994). 

15. S I Svinolupov. Second-order evolution equations with symmetries. Russian Mathematical Surveys 
40,241-242(1985). 

S. I. Svinolupov and V. V. Sokolov. Vector-matrix generalizations of classical integrable equations 
Theor. Math. Phys. 100, 214-218, (1994). 

16. V. Gerdjikov, A. Kyuldjiev, G. Marmo, G. Vilasi. Complexifications and Real Forms of Hamiltonian 
Structures. European J. Phys. 29B, 177-182 (2002). 

17. H. E. Nistazakis, D.J. Frantzeskakis, PG Kevrekidis, BA. Malomed, and R. Carretero-Gonzt'alez. 
Bright-Dark Soliton Complexes in Spinor Bose-Einstein Condensates. Phys. Rev. A 77, 033612 (2008). 

18. V. S. Gerdjikov. Algebraic and Analytic Aspects of A^-wave Type Equations. nlin.SI/0206014 ; Con- 
temporary Mathematics 301, 35-68 (2002). 

19. VS. Gerdjikov. The Zakharov-Shabat dressing method amd the representation theory of the semisim- 
ple Lie algebras. Phys. Lett. A, 126A, n. 3, 184-188, (1987). 

20. R. Ivanov. On the dressing method for the generalised ZakharovUShabat system. Nuclear Physics B 
694 [PM] 509U524 (2004). 

21. A. B. Shabat. Inverse-scattering problem for a system of differential equations. Functional Analysis 
and Its Applications, 9, 244-247 (1975). 

— An inverse scattering problem. Diff. Equations, 15 1299-1307 (1979). 



